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1 INTRODUCTION 



Star clusters are subjected to various disruptive effects, 
which prevent low mass star clusters (M c < 1O 4 M0) from 
surviving for a Hubble t ime. From o bservations thi s was 
first noted bv lOord Jl95gl) and later bv lWielenl Jl97ll) . who 
showed that there is a lack of open star clusters older than 
a few Gyr in the solar neighbourhood. Later, more quanti- 
tative results were obtained for the life time of star clusters 



ABSTRACT 

We investigate encounters between giant molecular clouds (GMCs) and star clusters. 
We propose a single expression for the energy gain of a cluster due to an encounter with 
a GMC, valid for all encounter distances and GMC properties. This relation is verified 
with ./V-body simulations of cluster-GMC encounters, where the GMC is represented 
by a moving analytical potential. Excellent agreement is found between the simulations 
and the analytical work for fractional energy gains of AE/\Eq\ < 10, where \Eo\ is the 
initial total cluster energy. The fractional mass loss from the cluster scales with the 
fractional energy gain as (AM/M ) = f(AE/\E n \), where / ~ 0.25. This is because 
a fraction 1 — / of the injected energy goes to the velocities of escaping stars, that 
are higher than the escape velocity. We therefore suggest that the disruption time of 
clusters, tdi s , is best defined as the time needed to bring the cluster mass to zero, 
instead of the time needed to inject the initial cluster energy. We derive an expression 
for idis based on the mass loss from the simulations, taking into account the effect of 
gravitational focusing by the GMC. Assuming spatially homogeneous distributions of 
clusters and GMCs with a relative velocity dispersion of <r cn , we find that clusters loose 
most of their mass in relatively close encounters with high relative velocities (~ 2 <7 cn ). 
The disruption time depends on the cluster mass (M c ) and half-mass radius (rh) as 
idis = 2.0 5 (M C /10 4 M Q ) (3.75pc/r h ) 3 Gyr, with S = 1 for the solar neighbourhood 
and S scales with the surface density of individual GMCs (£ n ) and the global GMC 
density (p n ) as S oc (E n p n ) _1 . Combined with the observed relation between rh and 
M c , i.e. rh oc M c , <dis depends on M c as idis oc Mc 1 ■ The index 7 is then defined 
as 7 = 1 — 3A. The observed shallow relation between cluster radius and mass (e.g. 
A ~ 0.1), makes the value of the index 7 — 0.7 similar to that found from observations 
and from simulations of clusters dissolving in tidal fields (7 ~ 0.62). The constant of 
2.0 Gyr, which is the disruption time of a 10 4 M cluster in the solar neighbourhood, is 
about a factor of 3.5 shorter than found from earlier simulations of clusters dissolving 
under the combined effect of galactic tidal field and stellar evolution. It is somewhat 
higher than the observationally determined value of 1.3 Gyr. It suggests, however, that 
the combined effect of tidal field and encounters with GMCs can explain the lack of 
old open clusters in the solar neighbourhood. GMC encounters can also explain the 
(very) short disruption time that was observed for star clusters in the central region of 
M51, since there p n is an order of magnitude higher than in the solar neighbourhood. 

(e.g. lHodgelll988l and Ide Griis fc Andera 2006 for clusters 
in th e LMC:lBoutlqukos fc Lamersll2003l for the SM C and 
M33; lLamers et al.ll2005l and iPiskunoy et al.ll20od for the 
solar neighbourhood and lGieles^tHnTl27)o"5r "for the central 
region of M51). From the theoretical side, the evolution and 
disruption of star clusters has been the subject of many 
studies using a variety of te chniques to solve the gravita- 
tiona l iV-body problem (e.g. |Ostriker 1 Spitzer fc Chevalierl 
Il972l IChernoff fc Weinberslll990t iGnedin fc Ostrikerlll997l : 
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Takahashi fc Porteeies Zward 120001 iBaumgardt fc Makiiiol 
20031). 

ILamers. Gieles fc Porteeies Zwartl l|2005) compared ob- 
servationally determined disruption times in four galaxies 
to the corresponding times following from the JV-body sim- 
ulatio ns of clusters in tidal fields of Baumgard t fc Makinol 
J2003ft and lPortegies Zwart et"aT] Jl998l.l2002l'l . Thev found 
that, based on the results of the simulations, the disrup- 
tion time (tdis) should depend on the cluster mass M c as 
t dis oc [M c /ln(M c )] 75 oc M c - 62 . The index (0.62), which 
is referred to as 7, was also determined from the observed 
age and mass distributions of different cluster populations 
jBoutloukos fc Lamersll2003h and the mean value from the 
observations is 7 = 0.62 ± 0.06. The value of tdis scales 
with the disruption time of a 10 4 Mq cluster as tdis = 
£4 (M c /10 4 Mq) 7 . Good agreement for £4 between simula- 
tions and observations was found only for the star clusters 
in the SMC. In the solar neighbou rhood, the observed i4 
for open clusters (£4 = 1.3 Gyr, see ILamers et al1l2005tl is 
a factor of 5 lower than expected from the disruption time 
due to the galactic tidal field, combined with a realistic stel- 
lar mass function and with stell ar evolution (ti = 6.9 Gyr, 
sec IBaumgardt fc Makinol I2003I) . An even larger disagree- 
ment of almost an order of magnitud e was found for td is of 
clusters in the central region of M51 llGieles et al.ll2005h . 

The fact that the observed disruption times are shorter 
could be a result of time-dependent external perturbations, 
such as spiral arm passages and encounters with giant molec- 
ular cl ouds (GMCs), which were ignored in the sim ula- 
tions oflBaumgardt fc Makinol (120031) . iTerlevichl Jl987l) and 
iTheunsI 11991 " showed that a single encounter with a mas- 
sive GMC_coulo^djsnnM 1 j^ Mq 
and iGieles. Athanassoula fc Portegies Zwartl i2006fi showed 
that in the solar neighbourhood passages of spiral arms con- 
tribute significantly to the disruption of open clusters. 

In this study we investigate the disruption of clusters 
by encounters with GMCs. We want to quantify whether 
GMCs can explain the difference between the observed val- 
ues of £4 and the ones following from the simulations without 
the time-dependent external effects. Moreover, would simu- 
lations including encounters with GMCs preserve the nu- 
merical value of the index 7 ~ 0.62, where the observations 
and previous simulations already agree? 

The paper is organised as follows: In § [5] we introduce 
the initial conditions for the clusters and GMCs and de- 
fine the parameters involved in a cluster-GMC encounter. 
In § [3| we give analytical calculations of the energy gain for 
a cluster due to various types of encounters. These analyt- 
ical formulae for the energy gain are verified with iV-body 
simulations and compared to the mass loss of the cluster in 
§ ^] The results of the iV-body simulations are used in § |S] 
to derive an expression for the cluster disruption time in en- 
vironments with different GMC densities. The conclusions 
and discussion are presented in § |S] 



2 INITIAL CONDITIONS 
2.1 GMC properties 

In this study we consider encounters between molecular 
clouds and star clusters. Since 90% of the total molecular 



gas mass in our Galaxy is in giant molecular cl ouds (GMCs) 
with mass M n > 10 4 M dSolomon et al.lll987l) and encoun- 
ters with low mas s clouds do not affect a star cluster much 
fe.g. IWielenlll985f) . we consider only the massive GMCs. In 
§|3wc show that the cluster energy gain due to an encounter 
with a GMC scales with the GMC mass squared (M n 2 ), sup- 
porting our assu mption to consider on ly the massive clouds 
(i.e. the GMCs). ISolomon et alJ lll987l) showed that there is 
a relation between the size and mass of galactic GMCs of 

M » = 540 (f:) 2 M °> « 

where M n is the mass of the GMC and R u is its radius 1 . The 
internal density profile is of the form p(r) oc r _1 . Eq. im- 
plies a constant mean surface density of E n = 170M Q pc -2 . 

We use th e spherically symmetric Plummer model 
JPlummedlllml) to describe the GMC, since it is mathe- 
matically convenient to use in analytical calculations (§|3J. 
The potential of this model is 

$(r) = -GM/^/r 2 +a 2 , (2) 

where G is the gravitational constant, M is the total mass 
and a is the Plummer radius. If we choose the Plummer 
radius of a GMC equal to half the radius from Eq. ^ e -g- 
a n = 0.5 R n , the resulting profile is very similar to what 
would follow from a density profile of the form p(r) oc r~ . 
In Fig.Qwe compare the descriptions for a Plummer model 
(Eq. |5J and constant surface density profile (Eq. for the 
potential in the top panel, the acceleration in the middle 
panel and the density profile in the bottom panel. The Plum- 
mer model has the advantage that the first derivative of the 
force is differentiable for each r, which is necessary for the 
A^-body simulations. 

2.2 Star cluster properties 

For the star cluster we also choose a Plummer model (Eq.[2]>. 
The number of particles is = 2048 with a total mass of 
M c = 1000 Mq and all stars the same mass ] resulting in a 
realistic mean stellar mass of 0.49 M Q (e.g. lKroupall200lD . 
The virial radius is set to 2 pc. For a Plummer model, the 
virial radius relates to the Plummer radius of the cluster 
as r v = 16a c /(37r). The cluster half-mass radius relates to 
a c as rh ~ 1.305 a c . Therefore, the values of a c and rh in 
physical units are 1.18 pc and 1.54 pc, respectively. The en- 
ergy gain of a cluster due to an encounter with a GMC 
will be compared to the total initial cluster energy, e.g. 
the sum of kinetic and potential energy. This is defined as 
Eo — —riGM c 2 /(2r^), with r\ ~ 0.4, depending on the clus- 
ter model. For a Plummer model E = (37r/64) GM c 2 /a c , 
hence r/ ~ 0.38. 

2.3 Encounter definitions 

Two particles with a relative velocity at infinity that is non- 
zero will experience a hyperbolic encounter. The motion of 
the reduced mass particle is schematically represented in 

1 Throughout this paper we use the subscript "n" (from nebula 
or nuage) to indicate the parameters of the GMC and "c" for the 
cluster. 
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Figure 1. Comparison between the radial profiles of the potential 
(top), acceleration (middle) and density (bottom) for a constant 
surface density GMC with M n = R n = 1, following from Eg. HI 
and p n oc r —1 (dashed line) and for a Plummer model with radius 
a n = 0.5 fin (full line). 
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Figure 2. Motion of the reduced particle in a hyperbolic en- 
counter. The impact parameter b and the velocity at infinity vq 
are denoted on the left side of the picture and the encounter pa- 
rameter p and maximum relative velocity Vniax are indicated at 
the point of closest approach. 

There exists an analytical expression for the energy gain 
of a cluster due to an head-on encounter and due to a dis- 
tant encounters. In the next section these relations will be 
extended to the full range of encounter parameters, i.e. in- 
cluding close encounters. The exact boundary between close 
and distant encounters is hard to define. Here we will dis- 
tinguish three types of encounters: (1) head-on encounters 
for which p — 0; (2) close encounters, where < p/a n < 5 
and (3) distant encounters, which is for p ^ San- 



Fig. |21 The impact parameter and the velocity at infinity 
are referred to as 6 and Vo, respectively. The distance of 
closest approach and the maximum velocity are called the 
encounter parameter p and Vmax, respectively. Conservation 
of angular momentum relates Vq and V ma x as 



V m 



vo- 



(3) 



Consider a GMC with a Plummer model according to Eq. 
and mass M n and a cluster with mass M c which is a point 
source. From conservation of energy and Eq. |3] it follows 
that b and p are related as 



2G(M n + M c ) 



1 + 



V l\/P 2 +«n 2 



(4) 



Note that b is now a mixed function of encounter and im- 
pact variables. In the next sections we will show that the 
encounter parameter p and Vmax determine the internal en- 
ergy gain. Since we want to perform iV-body simulations 
that result in an encounter with a certain desired p and 
Vmax, we use Eq. Eland Eq.0]to express p and Vmax in terms 
of b and Vo. In addition, in our simulations we will not start 
with the GMC at infinity, hence we will need to take into 
account the potential energy contribution at the beginning 
of the simulation. Since the final expressions for b and vo 
as a function of p, Vmax, M n and M c are quite complex, but 
follow from simple mathematical exercises, we do not give 
them here. 



3 CLUSTER ENERGY GAIN AND MASS LOSS 
DUE TO AN ENCOUNTER WITH A GMC 

3.1 Distant encounters 



Spitzer (1958;, hereafter S58) derived an expression for the 
energy gain of a cluster due to a distant encounter with 
a GMC. He assumed that the GMC is a point mass and 
that the encounter velocity is much faster than the internal 
velocity of the stars in the cluster, which is referred to as the 
impulsive approximation. The cluster energy gain can then 
be expressed in the encounter, GMC and cluster parameters 



AE= 1_(2GM» 

3 \P 2 Vmai 



M c r 2 . 



(5) 



Here r 2 is the mean-square position of the stars in the clus- 



iTheuns! <ll99ll) used combined hydro-dynamical and N- 
body calculations of a GMC-cluster encounter, to show that 
Eq. |S] is very accurate. He also notes the inconvenience of 
the ill defined r 2 parameter. In principle r 2 follows from the 
cluster profile, but this variable is very sensitive to variations 
of the positions of stars far from the cluster centre, hence it 
is not constant during the simulation because of scattering 

2 S58 makes no distinction between b, vo and p, Vmax- In ? HI we 
show that p and Vmax are the most convenient ones to use for 
calculating the energy gain. 
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of stars out of the cluster. Even in simulations covering a 
short time span, r 2 will increase due to stars that get at 
larger distances from the centre due to 2-body encounters 
in t he core. 

lAguilar fc White! (Il985l) have shown with TV-body cal- 
culations of encounters between two Plummer models with 
radii ai and a%, that the predictions of Eq. |S] hold when 
p ^ 5max(oi,a2). In §^]we examine with TV-body simula- 
tions in what regime the results of Eq. |S] are valid, and at 
what moment during the encounter the value for r 2 should 
be taken. 



3.2 Head-on encounters 

Gravitational focusing will deflect many cluster orbits which 
start with b > 5a n to encounters with p < 5a n , for which 
the relation of S58 (Eq. |SJ does not hold. To obtain an ex- 
pression for the energy gain for close encounters, we first 
explore the opposite of distant encounters, namely head-on 
encounters (i.e. p — 0). In those encounters the extent of the 
GMC has to be taken into account, since assuming a point 
mass GMC would severely overestimate the energy gain. The 
energy gain in close encounters should be an interpolation 
between the d istant encounter and head- on results. 

Following Binn ev fc Tremainel 1)198711 (Chap. 7.2, here- 
after BT87) we can derive the cluster energy gain for an 
head-on encounter between a cluster and a GMC, both de- 
scribed by a Plummer model. Let (R, z) be cylindrical coor- 
dinates, with the z-axis coinciding with the GMC trajectory. 
Due to symmetry, the stars in the cluster get a velocity in- 
crease only perpendicular to the trajectory of the GMC, of 



Av (R) = 



2GMnR 



V max (i? 2 +a n 2 ) : 



(6) 



where analogous to the case of the distant encounters we 
use Vmax for the relative velocity between the cluster and 
the GMC. Eq.EJcan be used to calculate the energy gain of 
the cluster: 



AE = n [Av(R)] 2 T,(R)RdR. 



(7) 



Here £(ii) is the surface density of the cluster obtained by 
projecting parallel to the line of motion of the GMC. Since 
the cluster is described by a Plummer model, £(i?) is 

E (,R) = ^! (fl 2 + ac 2 r 2 , (8) 

7T 

so the energy gain of the cluster due to an encounter between 
two Plummer models is 



Ai? = 



( 2GM n a c 

V Knax 



o (i? 2 + a n 2 ) 2 (i? 2 + a c 2 ) 



7 dR.(9) 



The result of the integration is a function of a c and a n and 
is equal to l/12a 4 when a n = a c — a. In that case, Eq. |S| 
reduces to the result found by BT87: 



A_Ebt = — 



1 / GM n 



3 V Vmax Oj 



M c 



(10) 



When On 7^ ac, the result of the integration is a function of 
a c and a n . Introducing a new variable x = a n/o c , the ratio 
of the GMC and cluster radii, Eq.[5]can be written as 



AE= 1 -' GM » 



3 V Vmax Oi. 



' M 12( X 2 ~ 1) ~ 6( X 2 + 1) ln( X 2 ) nu 

...... Mc ,[ ' 

where in the first factor we recognise the result for a c = a n 
(e.g. Eg, 1101 and the rest is a correction factor C(x) which is 
only a function of x- We can write a more general expression 
for the energy gain of head-on encounters 



Ai? = A£ BT C(x), 
with 

12( x 2 -l)-6(x 2 + l) ln(x 2 ) f nr / I 

1 for x = I 



(12) 



(13) 



The expression for C(x) is rather complicated, but for real- 
istic values of x (i- e - 2 < x < 10), C(x) can be approximated 
with high accuracy by C(x) — 2x~ 3 (see Fig.^J. 

With Eq. ll2l and Eq. ll3l we have a general expression for 
the energy gain of head-on GMC-cluster encounters, when 
both are described by Plummer models, valid for any ra- 
dius and mass of the cluster and the GMC. In § 14.31 we will 
confront Eq. I12l and Eo. ll3l with results from iV-body simu- 
lations. 



3.3 A general expression for the cluster energy 
gain 

The energy gain for a close encounter, i.e. an encounter for 
which the cluster moves through the cloud within a few a n 
from the GMC centre, should be a smooth interpolation be- 
tween the result for a head-on (Eq. I12H and distant encoun- 
ters (Eq.EJ. When replacing 4r 2 /p 4 in Eq.|51by C(x)/ac 2 , 
we get the result for head-on encounter fEo. 1121 . One way 
of connecting the two is to add a term in the denominator of 
Eq. |S] preventing the result to diverge to infinity at p = 0, 
so that it converges to Eq. 1121 is 



AE 



4r 2 



P 2 + 



^4r 2 a c 2 /C(x) 



_ ( GMn 
2 I V 

\ v max 



(14) 



If we use the approximation given in § 13.21 for C(x) (i-e. 



C(x) - 2x ) and the fact that Vr 2 ~ 1.8 rh, a c ~ r h /1.305 



and the real GMC radius R n relates to a n as R 
can rewrite Eq. 1141 as 

At: — t I — - I M c . (15) 



p 2 + y/ r h Rn 



(GMn 
\ Vmax 



In the next section, this relation will be compared to results 
of iV-body simulations for encounters with p ranging from 
to 10 a n . 



4 VALIDATION WITH TV-BODY 
SIMULATIONS 

4.1 Description of the code 

The TV-body calculations were carried out by the kira in- 
tegrator, which is part of the Starlab software envir onment 
jMcMillan fc HutHl996t iPorteeies Zwart et al-lboOlf) . Kira 
uses a fourth-order Hermite scheme and includes special 
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treatments of close two-body and multiple encounters of ar- 
bitrary complexity. T he special purpose GRAPE-6 systems 
(Mak ino et, alJ l2003:'l of the University of Amsterdam are 
used to accelerate the calculation of gravitational forces be- 
tween stars. 

The potential, acceleration and jerk of the stars due to 
the GMCs are derived from Eq. [5] and calculated for each 
star individually, based on its position and velocity with 
respect to the GMC. We assume the GMCs are moving 
potentials which do not interact. Indeed, we are interested 
in the star cluster and since we consider GMCs that are 
~ 100 — 1000 times more massive than the cluster, the grav- 
itational pull from the star cluster on the GMCs is negligible 
and will therefore not affect its orbit. Although GMCs have 
finite life times (~ 10 — 30 Myr), the encounter durations are 
sufficiently small (~ 1 Myr) to not be affected by this. 



4.2 Description of the runs 

We want to avoid sudden jumps in the contribution of the 
GMC to the cluster potential energy at the beginning and 
at the end of the simulations. In principle, our GMC should 
start at infinity, have an encounter with the star cluster 
and then continue until infinity. Since we do not want to 
evolve the cluster to long before the encounter, to restrict 
simulation time and to avoid too much internal dynamical 
evolution, we put the GMC initially at a distance of D — 
500 a c from the cluster and give it zero mass. While the 
GMC travels from D to 0.9 D, approaching the cluster, we 
increase its mass gradually from to M n . Similarly, after 
the encounter and while the GMC travels from 0.9 D to D 
we decrease its mass gradually from M n to 0. To give the 
cluster time to get back to equilibrium and to give stars 
with velocities higher than the escape velocity time to leave 
the cluster we evolve it further for a time equal to half that 
elapsed after the encounter. 

To test the different analytical expressions of § [3] for 
the energy gain of the cluster, we perform series of head-on 
(p — 0), distant (p > a n ) and close encounters with different 
GMC radii and encounter parameters. The initial velocity 
(vo) and impact parameter (b) of the GMC are calculated 
from the desired values of p and V max , using the expressions 
for gravitational focusing of § 12.31 The cluster properties are 
described in § 12.21 A summary of all runs is given in Table0 
The results will be discussed in § 14. 314.51 

All clusters are scaled to JV- body units, such that G = 
M — 1 and Eq — —0.25, following lHeggie fc MathieuHl986D . 
Then r v is the unit of length and follows from the scaling 
of the energy, since r v = GA/ c /(4|£'o |) = 1, where Eo is 
the total (potential and kinetic) initial energy of the cluster. 
We assume virial equilibrium at the start of the simulation, 
which implies r v (t = 0) = 1. 



4.3 Head-on encounters 

We perform iV-body simulations of head-on encounters be- 
tween a cluster and a GMC, both described by a Plummer 
model. GMCs with different Plummer radii a n were used 
to compare simulations to the analytical result of Eq. 1121 



The energy gain of the cluster 3 is always expressed in units 
of the initial cluster energy, i.e. the fractional energy gain 
(AE/\Eo\). Since we are here interested in the effect of the 
GMC radius on AE/\E \, we fix the mass of the GMC at a 
constant value. The results can then be compared with the 
predicted value from Eq. 1121 More specific, when we divide 
the energy gain for the different cases by the value that was 
predicted by BT98 (e.g. Eq.Hrjl. we should get the factor 
C(x) °f Eg. 1131 In § !4.5l we show that there is no good agree- 
ment between the simulations and the analytical predictions 
when AE/\E \ » 1. To keep AE/\E \ < 1 for all GMC 
radii, we choose M n = 10 4 Mq , which is somewhat lower 
than a realistic value for M n . We choose integer values for 
X from 1 to 10, where Eq. 1101 and Eo. ll2l should be valid for 
X = 1 and x / 1, respectively. The encounter velocity for 
all simulations is Knax = 20<r r ms, where <r rms = lkms -1 is 
the internal velocity dispersion of stars in the cluster. 

In Fig. 0] we show the results of A_B/|_Bo| from the N- 
body simulations for different x with diamond symbols, the 
value of BT87 with a filled circle (left) and the result of 
Eq. 1121 with the full line. The latter relation is also shown 
with a dashed line, but with an approximation for C(x) of 
the form C(x) = 2x~ 3 - This approximation is very close to 
the real C(x) for 2 ^ x ^ 9- Excellent agreement between 
the analytical work and the iV-body results is found in all 
cases. 

In Fig. [3]we show snapshots of a cluster during a head- 
on encounter with a GMC. In the top panels the line of 
sight is perpendicular to the trajectory of the GMC, whose 
direction is parallel to the a;-axis. The attractive force of 
the GMC parallel to its trajectory before and after the en- 
counter cancel out, resulting in displacement of stars only 
perpendicular to the line of motion of the GMC. In the bot- 
tom panels the line of sight is along the motion of the GMC, 
and it can be seen that the stars escape and form a disk like 
structure of unbound particles around the cluster. The time 
with respect to the encounter moment is indicated in each 
panel in the top row, in units of i onc , where t cnc = 2o n /V r ma x. 

4.4 Distant and close encounters 

For the simulations of distant encounters we choose M n = 
1O 5 M . The tidal forces of the GMC deform the cluster 
into a cigar like shape, as can be seen in Fig. |S] In Fig. |S] 
we show with the full lines the evolution of the internal en- 
ergy of the cluster for encounters with three different val- 
ues for p as a function of time, in units of the encounter 
time (tone = p/Knax). (Note that for head-on encounters 
t onc = 2a n /V r ma x.) We define time (t) such that the moment 
of closest approach is at t = 0. The value of Knax for these 
runs is 40<r rms and the GMC radius is equal to the cluster 
radius. The GMC radius is chosen small in this example in 
order not to have too large values of p when the ratio p/a n 
is large. Otherwise, the energy gain would be too small. The 
results of Eq. 0for these encounters are overplotted as the 

3 We start the cluster in the origin, resulting in a centre of mass 
velocity after the passage of the GMC. In the comparison with 
the analytical models, when we refer to cluster energy, we refer to 
the sum of the internal potential and internal kinetic energy, i.e. 
the energy after the centre of mass motion has been subtracted. 
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Figure 3. Snapshots of a star cluster undergoing a head-on encounter with a GMC of M n = 10 4 Mq with V max = 20<r rms . The motion 
of the GMC is along the x-axis the line of sight is perpendicular to (top panels) and along (bottom panels) the GMC trajectory. The 
arrows in the left lower corner of the left panels are parallel to the direction of motion of the GMC. The GMC is shown with grey shades 
based on the surface density of a GMC with a n = 5.8 a c . The time with respect to the moment of encounter is indicated in each panel 
of the top row, where t cnc = 2a n /V max . The snapshots are taken in the centre-of-mass frame of the cluster and the viewing distance is 
40 ac. 
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X = a„/a c 

Figure 4. Fractional energy gain (AE /\Eo\) for a star cluster 
due to head-on encounters with GMCs with M n = 10 4 Mq for 
different \ = a n / a c and with Vmax = 20 0" rms • The analytical 
approximation of BT87 (Eq. I1UI for head-on encounters between 
identical Plummer models is indicated with a filled circle (left). 
The iV-body results are shown with diamonds. The analytical 
prediction of Eg. I12l with the correction factor C(\) for different 
cloud radii is shown with the full line. The predicted energy gain 
with an approximation for C(\) (i.e. C(x) = %X~ 3 ) i s shown 
with the dashed line. 



dashed lines. The value of r 2 changes during the simulations 
due to internal dynamical effects and because of deforma- 
tion of the cluster by the GMC. This last effect is clearly 
visible in Figs.|3and|S| Since the result of Eq.[K]is very sen- 
sitive to the value of r 2 , we performed tests to find at which 
time r 2 should be taken to get the best agreement between 
Eq. and the simulations. If we take this value too early, 
Eq. will underestimate the energy gain compared to the 



result of the simulation. This because dynamical evolution 
will increase r 2 . When using a value r 2 which is determined 
during the encounter (i.e. t ~ 0), the energy gain from Eq.0 
overestimates the true value. We found that best agreement 
was found when using r 2 at t — —2t euc . From Fig- El it can 
be seen that this is moment where the energy of the clus- 
ter starts to increase, whereas the cluster energy is almost 
constant before that moment. 

The variations of the relative velocity between the clus- 
ter and the GMC are smaller than the variations in r 2 , and 
for that reason the result of Eq.Q)]is less sensitive to the time 
that we choose the velocity. Since most of the energy is in- 
jected between — t cnc and +t cnc (Fig.|SJl, the best agreement 
is obtained when the relative velocity at t = 0, e.g. Vmax, is 
used. 

To test the validity of the interpolation between the 
distant and the head-on encounters (Ea. I14jl . we performed 
JV-body simulations of various encounters with impact pa- 
rameter p ranging from to 10 times the GMC radius o n . 
The encounter velocity was chosen 40<r rms , such that the 
fractional energy gain is always lower than 1 (Ea. II 2p> . The 
impact parameter b and initial velocity Vq of the GMC where 
calculated using Eq. |1] The GMC mass is 1O 5 M and we 
used the observed mass-radius relation for GMCs (Eq^l , re- 
sulting in x = 5.8 and therefore C(x) = 0.01 (Eq. In 
Fig. Q we show the result of AE/\Eo\ from the iV-body sim- 
ulations (diamonds) and the analytical prediction of Eq. 1141 
(full line). The predictions for the two extreme values of p, 
i.e. for distant and for head-on encounters are shown with 
the dashed line and the filled circle (left), respectively. The 
analytical expression Eq. 1141 agrees very well with the re- 
sults from the JV-body simulations. It can be seen that the 
prediction of S58 is accurate for p ^ 5 a n , in a greement with 
what was found by lAguilar fc White! (119851) . Note that in 
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Table 1. Summary of the JV-body simulations. All runs mentioned are done 4 times and average results are used. (1) GMC mass in 
units of M c = 10 3 Mq; (2) ratio of the GMC radius and cluster radius (x = a n /a c , where a c =1.18 pc); (3) impact parameter in units 
of GMC radius; (4) impact velocity in units of internal cluster velocity dispersion (cr rms ); (5) and (6) are the encounter distance and 
velocity, following from columns (3), (4) and Eq. [I]and Eq. [3] (7) Figure numbers where the results of the runs are shown. 
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Fig. 0we have used scale length to normalise p. This 

is possible since GMCs have finite sizes. If we were dealing 
with objects of infinitesiaml size, as e.g. black holes, this 
scaling would not be possible. Thus, Fig. can not be ap- 
plied to those cases. For large p (p » a n ) the energy gain 
only depends on p. This is demonstrated by the good agree- 
ment with the point source approximation of S58, since for 
these cases a n = 0. 

To see if the agreement between our predicted energy 
gain from Ea. ll4l and iV-body simulations holds for lower ve- 
locities, we performed a similar series ./V-body simulations 
as shown in Fig.Qfor Vm ax = 10er rms and 20<r rms . The values 
of AE/\Eq\ for these three different values of V m ax and dif- 
ferent values of p is shown in Fig. [8] For the lower velocities 
we again find very good agreement. However, when the frac- 
tional energy gain is larger than 1, the iV-body simulations 
show a smaller energy gain than the predictions. The dif- 



ference is la rger for larger ener gy gains. Similar results were 
obtained bv lGieles et al.l (1200611 for one-dimensional shocks 
due to spiral arms. This is probably because the impulsive 
approximation does not hold anymore. When the GMC is 
very massive, or when the relative velocity is low, the stars 
in the cluster are displaced too much before the encounter. 
The symmetry assumption of the encounter of the impulse 
approximation then does not hold anymore. 



The energy gain from the encounter is absorbed mainly 
by the stars far from the cluster centre, which can gain ve- 
locities much higher than the escape velocity. A bound core 
will remain, even when AE/\Eq\ >> 1. Both topics will be 
discussed in more detail in § 14.51 
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Figure 6. Snapshot of a star cluster during a distant encounter with a GMC with M n = 100 M c , p = 3a n and Vmax = 15<r rms in 
the centre of mass frame of the cluster. The motion of the GMC is along the rc-axis, starting at b = y = 20 a c , and the line of sight 
is perpendicular to (top panels) and along (bottom panels) the GMC trajectory. The arrows in the left lower corner of the left panels 
are parallel to the direction of motion of the GMC. The GMC is shown with grey shades based on the surface density of a GMC with 
a n = 5.8a c . The time with respect to the moment of encounter is indicated in each panel of the top row, where t enc = p/V m ax- The 
snapshots are taken in the centre-of-mass frame of the cluster and the viewing distance is 60 a c ■ 



TV- body 
Spitzer (1958) 




p - 20a n 




\ Spitzer (1958) 



<^> N-body 
• AE m * C(x) 

Eq. 15 




Figure 5. The evolution of the energy gain of a cluster dur- 
ing distant cloud encounters (full lines) with M n = 2 X 10 5 Mq, 
Vmax = 30 <r rms and three different impact parameter p. The 
dashed lines indicate the final energy gains predicted by S58 
(Eq.|SJ. 



4.5 Relation between predicted energy gain and 
mass loss 

In the previous section we found good agreement between 
the TV-body simulations and our analytical formula for 
Ai?/|i?o| of a cluster due to encounters with GMCs with 
different p and Vmax as long as AE/\Eo\ < 1 (see Eg. 1141 . In 
Fig. |H1 we show that deviations occur whenever the energy 
gain is larger than the initial cluster energy. We did a few 
more simulations of head-on encounters with higher GMC 
masses and compare all the predicted energy gains to the re- 
sults from the simulations. The results are shown in Fig. |2l 
A one-to-one relation (dashed line) holds until AE/\Eq\ ~ 5. 



Figure 7. Fractional energy gain of a cluster as a function of the 
impact parameter p and for Vmax = 40 o- IulB . The results of the N- 
body simulations are shown with diamonds. The analytical result 
for head-on encounters, including the correction for the radius of 
the GMC is shown with a filled circle (left). The result for distant 
encounters from S58 is shown with the dashed line. The result 
from the interpolation between the two analytical predictions (e.g. 
Eq.l51and ll2l . i.e. Eg. 1141 is shown with the full line. 

A simple functional (y — x[l — exp(— 12/:r)]), which starts 
to deviate from a one to one relation around y ~ 5 and lev- 
els off at y = 12 is also shown in Fig. [5] and describes the 
results well. 

To derive a realistic disruption time of clusters due to 
encounters with GMCs, we need to know how the energy 
gain relates to the number of stars, or mass, that is lost from 
the cluster in the different encounters. In Fig. 1101 we show 
the fractional mass loss (AM/Mq) as a function of AE/\Eo\. 
The dashed line indicates a one-to-one relation and the full 
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Figure 8. Fractional energy gain of a cluster as a function of 
the impact parameter p in units of a n and for different values 
of Vmax- The results of the TV-body simulations are shown with 
diamonds, where the size of the diamonds corresponds to the 
encounter velocity. The analytical results from Eg. 1141 are shown 
with the full lines. 




1CT 4 10~ 2 10° 10 2 

^-^ANALYTIc/ ^0' 



Figure 9. Relation between the fractional energy gain from the 
simulations and the analytical prediction of Eq. 1141 (circles). A 
one-to-one relation is overplotted as the dashed line. A simple 
functional form, describing a linear function that saturates at a 
value of 12, fits the result at the high energy part well. 



some encounters with A E > 10i?o, called overkill encoun- 
ters llWielenlll985U l988). can completely disrupt a cluster. 
We will further discuss this topic in § 15.31 

Since we consider isolated clusters, one might argue that 
there can be stars that have remained bound, but that would 
have been pushed over the tidal radius if the cluster was lo- 
cated in a galaxy. Assuming a disk galaxy with a fiat rotation 
curve, with rotational velocity of Vo = 220 km s -1 , the tidal 
radius of a cluster depends then on its mass and galacto- 
centric distance (7?) as r t = (GM C /2V 2 ) 1/3 R 2/3 . The esti- 
mated values of r t for our cluster at R = [0.5,1,2] x Ro, 
where Rq — 8.5 kpc, i.e. the solar neighbourhood, are 
r t ~ [13.8, 21.8, 34.6] pc, respectively. We counted for our 
simulations the number of extra stars (AN /No) which af- 
ter the encounter beyond these three radii. There is a lin- 
ear relation between AN /No and AM /Mo, where the ratio 
of the two quantities is [1.3, 1.1, 1.0] for the three different 
tidal limits, respectively. So the effective value of / is some- 
what higher when we take into account the tidal limit of 
the cluster. For the solar neighbourhood a good estimate 
would be / = 0.22 x 1.1 ~ 0.25. For clusters in the centres 
of galaxies the tidal radius is small and / can be a factor 
~ 1.5 higher. From now on we will assume a linear relation 
AE/\E \ = / AM /Mo, with ~ 0.25, which takes into ac- 
count stars that have become unbound in isolation and/or 
are pushed over the tidal limit. This will be used in §|3]to 
derive the mass loss of clusters in environment with GMCs 
and the resulting cluster disruption time. 

Hitherto, the cluster disruption time has always been 
defined as the time needed to bring the cluster energy 
to zero by pe ri odical ly injecting energy with shocks (e.g. 
lOstriker et all (Il972ft for the disru ption time by disk 
shocks and lSpitzerfcCnevalierl dl973l) for the disruption by 
GMCs). We have shown that this overestimates the disrup- 
tion time by a factor of 1/f ~ 4. A more accurate disruption 
time, taking into account this factor / is derived in the fol- 
lowing section. 



5 CLUSTER DISRUPTION TIME DUE TO 
ENCOUNTERS WITH GMCS 

With the description of the mass loss from a cluster due 
to different types of encounters with GMCs of §§ 13 141 we can 
derive the disruption time of star clusters in an environment 
with GMCs, such as the solar neighbourhood. 



line is a fit to the results. The grey areas indicate the regions 
where AM /Mo cannot be defined, since the cluster can not 
loose less then one star (i.e. Mo/2048, lower region), and it 
can not loose more than the total initial cluster mass (i.e. 
Mo, upper region). The fit shows an almost linear relation, 
where AM /Mo is a fraction / = 0.22 of AE/\Eo\. Given the 
good agreement between the predicted energy gain and the 
one from the simulations, we can conclude that a fraction 
/ = 0.22 of the input energy is used to unbind stars. This 
implies that a fraction (1 — /) of AE/\Eo\ went into the 
velocity of the escaping stars, so that these have escaped 
with velocities much higher than the local escape velocity. 
Combined with the saturation of Fig. |5J this shows that 



5.1 Clusters in an environment with GMCs 

When considering a distribution of GMCs and clusters with 
a relative velocity dispersion a cn = \/a'l + o'£ , the proba- 
bility of a relative velocity in the range vo,vo + dvo can be 
expressed as (BT87, Eq. 7-60) 



dP = 



4nvo 



(47TCr cn 2 )3/2 



exp 



4<7 cn 



dv . 



(16) 



Assuming a spatially homogeneous distribution of GMCs, 
the rate of encounters in the range b, b + db is then 



C 



ridouds Vo 2n b db dP 



(17) 



10 M. Gieles et al. 



10.0000 



1 .0000 



5g 0.1000 

< 

0.0100 



0.0010 



0.000 

0.0001 



r - - 


y - * 

y = 0.22»x - 96 






/ 






/ 0/ 










y ' , A,,,A I I I I 



0.0010 



0.0100 
AE, 



000 
n-body/'-^'o 



.0000 10.0000 



Figure 10. Relation between the fractional mass loss (AM/Mq) 
and the fractional energy gain (AE/\Eq\) for the simulations of 
§ |3j (circles) . The dashed line shows a one-to-one relation and the 
full line is a fit to the simulations. The grey areas indicate the 
regime where AM/Mq can not be defined. 



^clouds 87T 2 b db 

(47RT cn 2 )3/2 



exp 



4<7 cn 



duo, 



(18) 



where n c i ou d s is the number density of clouds. The energy of 
the cluster increases at a rate 

P OO I' b ma x 

E= / C(b,v ) AE(b,v )dbdv , (19) 

where AE(b, Vo) follows from Eg. 1141 after converting b and 
wo to p and V m ax due to gravitational focusing. This integral 
can be calculated numerically. 

In Fig. 1111 we illustrate this using the parameters for 
clusters and GMCs used in the iV-body simulations of § ^ 
e.g. M n = 10 5 M Q , a n =6.8pc and vV ~ 1.8 pc. We assume 
that young clusters and GMCs have the same velocity dis- 
persi on, i.e. <j n = a c — 7kms _1 l)St,arkll984l : IPiskunov et alJ 
2006). The velocity dispersion of the motion of GMCs with 
respect to clusters, is then a cn — \/a'l + a'£ ~ lOkms -1 . 
In the left panel of Fig. 1111 we show the impact rate C of 
Eg. 1181 in units of ra c i ou d s as a function of b and vo, in units 
of a n and a cn , respectively. In the middle panel we show 
AE/\Eo\ for individual encounters as a function of b/a n and 
"Uo/ccn, where we have converted b and vo to p and vo first, 
and then substituted these values in the equation for en- 
ergy gain (Eo. 1141 . The rate of energy gain as a function 
of b/a u and wo/o"cn is the product of the encounter rate 
times the energy gain per encounter and is shown in the 
right panel. The isocontours indicate the regions where 10, 
50 and 75% of the total energy per unit time is produced. 
This total energy was not determined from this graph, but 
from a numeric integration of a similar figure, but with an 
area that is 100 times bigger. Most of the energy is injected 
by fast (vo/a cn ~ 2) and close (b/a B ~ 1) encounters. The 
encounter velocities (Vmax) are even higher than 2a cn due to 
gravitational focusing and this means that the impulse ap- 
proximation is valid for these encounters. In § 14.41 we found 
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Figure 12. The correction factor g for close encounters (see text 
for details) as a function of GMC mass. The full line shows the 
result when gravitational focusing is ignored, and the other lines 
show the results for 3 different values of cr cn when gravitational 
focusing is considered. The predicted value of g = 3 of BT87 is 
indicated with the dashed horizontal line. 



that the distant encounter regime of S58 holds for p > 5 a n . 
The peak in the right panel of Fig. 1111 around b ~ a n il- 
lustrates that the approximations of S58 for a point source 
GMC are not applicable for the encounters that inject the 
majority of the energy. In that regime it is important to use 
the more accurate description for the energy gain found in 
SETSlrEa.lTH. 

BT87 solved Eq.^Jby substituting the energy gain for 
distant encounters of S58. They then integrate b from R n to 
infinity and correct for the missing close encounters with a 
factor g of order unity. The advantage of this solution is that 
the result has R n 2 in the denominator, which can be com- 
bined with M n and n c i ou d s to express the final result in the 
observable molecular gas density (p n ) and surface density 
of individual GMCs (S n ). Our results allow us to quantify 
the parameter g of BT87 by integrating Eo. 1191 numerically 
from to infinity and from R n to infinity. When gravita- 
tional focusing is ignored, the value of g depends only on 
M n , since in that case p — b and V ma x = vo- Including grav- 
itational focusing makes g dependent on a cn as well, since 
then encounters with initial b > R u can result in p < R n . 
Intuitively one thinks that the inclusion of gravitational fo- 
cusing will increase the value of g, since there are more 
encounters with p < R n to correct for. However, the en- 
counter velocity V max is always higher than Vo, which will 
make the energy gain of the cluster for these encounters 
lower. In Fig. 1121 we show the resulting g for various val- 
ues of Mn ignoring gravitational focusing (full line) and for 
three values of a cn . The horizontal dashed line indicates the 
estimated value of BT87 for g. For the solar neighbourhood 
(e.g. M n ~ 10 5 Mo,(j C n — lOkms -1 ) our results show that 
g ~ 2.5, slightly lower than predicted by BT87. 

Using the relation between energy gain and mass loss 
of § 14. 51 we can define the mass loss rate of the cluster as 



M c 



j E 



(20) 



where / = 0.25 is a dimensionless constant that relates frac- 
tional energy gain to fractional mass loss (see Fig. E3 and 
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Figure 11. Left: The encounter rate (C, Eg, 1181 . relative to rt c i ouds , for different values of b and no- Middle: Fractional energy gain 
fEq, 1141 for encounters with different b and vg. Right: Rate of energy increase, resulting from the multiplication of C and AE. 



§ 14.51 . The initial cluster energy is E = -yGM 2 /(2r h ) 
wher e rj ~ 0.4, depending on the cluster profile iSpitzerl 
Il987l) . 

We can derive a relation for the mass loss, using the ex- 
pression for E of BT87 (their Eq. 7-66) , where we substitute 
E„ = M n /(irRn 2 ), pn = Mn n c ioud 8 and M — f E 



M c = 



-^—-GPnPn,^ 



(21) 



where the ratio (r 2 /r n ) depends on the cluster model. For 
the Plummer cluster that we used (r 2 /Vh) — 4 and for 
King profiles with dimensionless potential depth values of 
Wo = [3,5,7,9] it follows that (r 2 /Vh) ~ [1.5,2,3.5,3.7]. 
From Eq. Hj] we can derive the disruption time of clusters, 
by defining 



tdis = M c /M c 



3a cn ri 



Mc 
r h 3 ' 



(22) 



For the constants in Eq. 1221 we have to make several as- 
sumption. First, we assume / = 0.25, as was found in 
§ 14.51 For (Ton we take lOkms -1 , b ased on observations 
of GMCs in the s olar neighbourhood iStark fc Brandl[l989l) 
and open clusters dPiskunov et al 120061) . For the cluster pro- 
file we assume a King profile with Wo = 3, corresponding 
to c e log (r t /r c ) — 0.7, whi ch is the m ean concentration 
of the galactic open clusters Il962l) . King profiles are 

more realistic due to the tidal truncation, which the Plum- 
mer model does not have. For the Wo = 3 profile considered 
here, (rg/r 2 ) ~ 0.67 and rj = 0.42 (BT87). 



5.2 Dependence on the relation between cluster 
mass and radius 

The solution for tdis of Eq. 1221 is very sensitive to the choice 
of rh. Observations of (young) extra-galactic star clusters 
show that the projected half-light radius (r e «) of the major- 
ity o f star clusters is confined to a na rrow range of 2-4 pc 
(e.g. iBastian et al] 120051: iLarsenl [20041 and a very shallow 
relation between the radius and the mass. A factor of two 
difference in radius and nearly equal masses still results in 



a factor of eight difference in tdis- This implies that clus- 
ters that are a few times larger than the mean radius are 
very unlikely to survive and cluster that are a few times 
smaller will be almost insensitive to encounters by GMCs. 
We choose the use the recently ob served relation between 
radius and mass of ILarsenl (|2004), who fits a function of 
the form r c s — AAI C X to a sample of young star clusters in 
nearby spiral galaxies and finds that A — 1.12 ± 0.35 and 
A = 0.10 ± 0.03. There is a large scatter around the fit rela- 
tion, which, combined with the high sensitivity of tdis on rh 
(Eq. 1221 . makes it hard to r elate mass to radius. However, 
since the relation of ILarsenl i2004l) reflects the median rela- 
tion between r e ff and M c , the conversion will probably hold 
for a large cluster population. To convert r g to rh we use 
the relation rh = 4/3 r c ffjjSpjt zerj jl98 7l) and with that we 
can rewrite the result of ILarsenl a2004f) as 



r h = 3.75 



Mc 



(23) 



. 10 4 M ( 

with A ~ 0.1. Substituting Eo. 1231 in Eq. 1221 and inserting all 
constants, yields an expression for the disruption time that 
is independent of rh and depends on the assumed index that 
relates rh and M c 



t dis = 2.0 



5.1 M, 



o 2 pc" 



En /On 



M c 



10 4 M Q 



Gyr, 



(24) 



with 7 = 1 — 3A. We normalised the product E n p n to 
the value for the solar neigh bourhood (S n = 17OM0pc -2 
and p n = 0.03M^pc~ 3 . ISolomon et al.lll987l) . When taking 
the 1 sigma upper bound value of A ~ 0.13 from ILarsenl 
(2004), the index 1 — 3A = 0.61 agrees v ery well with 
the v alue of 0.62 found from observations jLamers et alJ 
2005) and from TV-body simulations of clusters in the so- 
lar neighbourhood, dissolving under the combined effect of 
a tidal field, a realistic mass function and stellar evolution 
jBaumgardt fc Makinol 120031: | Taka hashi fc Portegies Zward 
l200ol) ~These simulations adopted a slightly different rela- 
tion between rh and M c since they assumed that at a given 
galactocentric distance the cluster size scales with the tidal 
radius (rh oc Mc 1 ' 3 ). The multiplicative constant of Eo. 1241 
is a factor of about 7 higher than what was found from ear- 
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lier studies (e.g. S58: ISpitzer fc Chevalierlll973l: BT87). The 
difference can largely be explained by our factor / ~ 0.25 
from § 14.51 that relates energy gain to mass loss. The re- 
maining difference is caused by the lower value we found for 
g and a slightly different value for the product E n p n . 



5.3 Disruption by single encounters 

In the derivation of tdis in the previous sections we integrated 
over the full range of b and v o in Eq. 1191 to get an expres- 
sion for the rate of energy gain of a cluster due to GMC 
encounters. Doing so, we implicitely assume that disruption 
is caused by a large number of encounters. However, if the 
cluster is completely disrupted by the first encounter with 
a GMC, Eg. 1221 will overestimate tdis- Fig. 1101 shows that 
a single enco unte r can indeed di srupt the cluster (see also 
IWielenlll985l and lTerlevicblll987lh If cluster disruption was 
always caused by single encounters, then Eqs. I22&I24I would 
overestimate the disruption time in all cases. The most accu- 
rate approach is to include the saturation function of Fig. [5] 
in the relation between energy gain and mass loss CEq.CT. 
i.e. to exclude encounters with AE > 10i?o. This results 
in AM/M = 1 for AE/\E \ > 10. However, this would not 
have allowed us to come to the convenient linear dependence 
of tdis on p n and £ mo i in Eq. |2"21 

Alternatively, we here estimate the typical time scale for 
a cluster to be di s rupte d by one single encounter (t d 'is Slc ). 
Following IWielenl £l985), we can calculate the critical en- 
counter distance (p cr it) for which one encounter will result 
in total disruption of the cluster. Based on Eq. 1151 on the 
expression for the cluster energy |i?o| = ??GM c 2 /(2rh) and 
on the relation between AE/\E \ and AM /Mo from § 1131 
we can derive p cr it by solving for p in A_E/|i5o| = 1//, such 
that 



Pcrit 



fSGf: 



V 3r? 



1/2 



V 



Mc 



-1/2 



- K. 



(25) 



In here K = \/ r n R n 3 , which is the term in the denominator 
of Eq. I15l that makes AE c onverge to the correct value for 
head-on encounters (p = 0) . IWielenl lll985l) did not include 
this term, since he assumed point-mass GMCs (Eq.0l. (Note 
that Eo. 1151 is equal to Eq. Qjjwhen K = 0). The number 
of encounters with encounter parameter p depends on the 
GMC number density (n c i ouds ) as 



N — Tvp n c i ou d s V At. 



(26) 



We can solve for At by setting N — 1 and p = p cr it, which 
is then the time it takes for a disruptive encounter to occur 



f lnsle = - . (27) 

dlS 7rpci-it 2 ^clouds V 

In her e gravitational focusing is ignored. Following IWielenl 
( 1985) and assuming point-mass GMCs, e.g. K = 0, we find 



an expression of e of the form 



f 

'■vi 



( St) 



' h 



\8Gfn 2 r 2 



1/2 



1/2 



(28) 



where we have used p n = M n n c i ou ds again. At first sight 
this relation for the disruption time looks similar to what 
we derived in the previous section (Eq. 12211 . The main dif- 
ference is the scaling with the cluster density (M c /r n 3 ). In 
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Figure 13. Comparison between tdis of Eq. 1221 (full line) and the 
single encou nter dis r uption time (tj'g 86 ). dotted line shows 
the result of Wielen 1 1985) for point-mass GMCs and the dashed 
and dotted-dashed lines show the results for extended GMCs with 
M B = 10 5 Mq and 10 6 M©, respectively. 



Fig. EDwe show a comparison between t d f B sle of Eq. HO (dot- 
ted line) and tdis of Eq.|22](full line) for different M c . Here 
we used a constant r n =3.75 pc and the same values for /, 77 
and (r 2 /r^) as in §0 The value of td^ glc is smaller than the 
result of Eq.^Jfor clusters with M c > 5 x 10 3 M Q . This im- 
plies that for massive clusters the life time is limited by the 
first disruptive encounter, rather than the successive heating 
by multiple encounters. For cluster of lower mass, the time 
it takes to disrupt by encounters with p > p cr it is shorter 



than t 



single 
dis 



This mass dependence is counter intuitive at 



3 in Eq.l2"71 based on M n = 10 5 M 
10 6 M (dotted-dashed line) and 



first, but is explained by the difference in scaling with M c 
(Eqs. HUMPS ■ Note that the crossing point (Fig. [T3J does 
not depend on p n , since Eqs. I22M28I depend in the same 
way on p n . 

In Fig. ll3l we also plot the result of td^ sle by substituting 
Eq. |31 with K = \/r h R, 
(dashed line) and M n = ru iviq 
Eq. Q Here we use the same value for p n , and then as- 
sume that all mass is in clouds with M n — 10 Mq and 
10 6 M , respectively, to solve for n c i ou d s (Eq. 1271 . For low 
mass clusters, these solutions for t B d ^ B slc agre with the re- 
sult of IWielenl jl985lh since there p cr i t > fewi? n , where the 
assumption of point-mass GMCs is valid (see Fig |7J . For 
clusters with 10 2 < M C /M Q < 10 4 , the solution for t d ^ gle is 
higher than that of the point-mass assumption. This is be- 
cause p cr it is smaller when extended GMCs are considered. 
In close encounters the extent of the GMC, e.g. K, is more 
important (Eq. I14H . There is no solution for clusters with 
M c > 5 - 10 x 1O 4 M , meaning that it is not possible to 
disrupt these clusters by a single encounters with a GMC 
with masses up to few times 10 6 Mq . 

The fact that the solutions for td^ Ele for realistic (ex- 
tended) GMCs are always larger than tdi s from Eq. [^justi- 
fies our assumptions of S I5.1I However, Fig. lK-il shows that for 
a cluster with M c = 10 3-4 Mq the two solutions are close, 
meaning that clusters are destroyed by just a few encoun- 
ters. This means that individual clusters will have strongly 
varying lifetimes. Our results of § [^1 will, however, be valid 
as a statistical mean. 
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6 CONCLUSIONS AND DISCUSSION 



We have derive an expression for the cluster disruption time 
(tdis) due to encounters with GMCs, based on the results 
of iV-body simulations of different types of encounters be- 
tween a cluster and a GMC. The expression for t dis (Eq.HU 
allows us to estimate the disruptive effect of GMCs in dif- 
ferent environments and to compare it to other effects such 
as the galactic tidal field, stellar evolution and spiral arm 
shocks. From Eg. 1241 we see that the disruption time for a 
10 4 Mq cluster in the solar neighbourhood is tn = 2.0 Gyr. 
This is a factor of 3.5 sho rter than what was found by 
iBaumgardt fe Making! <|2003T) for the combined effect of a re- 
alistic mass function, stellar evolution and the galactic tidal 
field (ti = 6.9 Gyr). Note that the index 7 is comparable, 
if we take the upper bound value of A = 0.10 ± 0.03, since 
Baumgardt fe Making! J2003F) found t dis oc [M c / ln(M c )] - 75 . 
Lamers et all (120051) showed that [M c /ln(M c )] 75 oc M C 0S2 . 
The same mass dependence was also found from the age 
an d mass distribution of c l uster samples in different galaxies 
bv lBoutlgukgs fe Lamersl i2003l) and from th e age distribu- 
tionof clusters in the solar neighbourhood bv lLamers et alJ 
(2005). For this last sample they found a shorter scaling 
value of ti = 1.3 Gyr. Eg. 1241 combined with the observed 
relation between and M c (Eq. 1231 suggests that the com- 
bined mass loss by GMC encounters and evaporation in a 
tidal field preserves the index 7 ~ 0.62. A similar scaling 
between tdis and M c was found for the disruption time by 
spiral arm passages. The combination of all these effects and 
a comparison to the observed age distribution of open clus- 
ters in t he solar neighbou r hood will be discussed in more 
detail bv lLamers fe Gieles! (I20061) . 

The result of Eq. 1241 can easily be applied to environ- 
ments/galaxies with different molecular gas densities (E mo i). 
The mid-plane density (p n ) of GMCs scales with E mo i as 
p n ~ E mo i/2/i, where h is the vertical scale length of molec- 
ular gas, which is a round 100 pc in the solar neighbourhood 
iDame et al.lll987l) . Assuming that h is the same in M51, 
and therefore p n = E mo i/200pc, we can predict £4 due to 
GMCs encounters once we know E mo i in M 51 and com- 
pare t his to the value of £4 f or clu sters in M51 dGieles et alJ 
l2005h . lGarcia-Burillo et al.l Jl993h find E mo i = 90 M Q pc" 2 
for the central region of M5 1, which is 14 times h igher than 
in the solar neighbourhood llSolomon et al.l [*1987). This im- 
plies ti ~ 2100 Myr/14 = 150 Myr. Note that this value is an 
upper limit, since we have assumed that the surface density 
of individual GMCs (E n ) is independent of E mo i, while this 
could be higher in environments with higher E mo i. In addi- 
tion, the values of / and g in Eq. I22l will probably be higher 
in the centre of M51 (see § !4.5l and §|HJ, which will make t d i s 
shorter. Also, we have assumed the same density profile as 
for the Galactic open cluster, e.g. Wo = 3. A density profile 
of Wo ~ 7, as found for the Galactic globular clusters, would 
reduce t d i s by a factor of two (Eq. . The resulting value of 
£4, however, is very close to th e value of £4 = 100 — 200 Myr, 
determined observationally bv lGieles et alJ i2005T) from the 
age and mass distribution of star clusters in M51. Based on 
this work we think that we can explain the short disruption 
times of clusters in the centre of M51, if GMC encounters 
are the dominant disruption effect. 
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